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Abstract 

o 

We investigate the confinement properties of solutions of the aggrega- 
tion equation with repulsive-attractive potentials. We show that solutions 
y—i remain compactly supported in a large fixed ball depending on the initial 

^ data and the potential. The arguments apply to the functional setting of 

. • probability measures with mildly singular repulsive-attractive potentials 

. ^ and to the functional setting of smooth solutions with a potential being 

X the sum of the Newtonian repulsion at the origin and a smooth suitably 

growing at infinity attractive potential. 



1 Introduction 



In this paper, we want to address "confinement" properties of solutions to the 
nonlocal interaction equation 



Pt = V-(p(Viy*p)), xeR^^>o, (1) 
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with compactly supported initial data p(0, ■) = po in a functional space to be 
specified. These nonlocal equations appear in many instances of mathematical 
biology [311 EH EHl ESI [23 E2], mathematical physics HH [HI ESI EZ] , and ma- 
terial science [HEllESllSHlElEESlIISlEn]. They are minimal models for the 
interaction of particles/agents through pairwise potentials. 

We say that the nonlocal equation ([T]) satisfies a confinement property in 
certain functional setting if every solution p{t, ■) to ([T]) in that setting with 
compactly supported initial data po is compactly supported for all times and its 
support lies in a fixed ball whose radius only depends on po and W. 

In most of the mentioned applications, particles/agents repel to each other in 
a short length scale while there is an overall attraction in larger length scales. 
Therefore, we can typically concentrate on repulsive-attractive potentials W 
as in [26liIl|22lllOl[25l[3Hl[33l|2Ol|2ll0 These potentials lead to a rich 
ensemble of compactly supported steady states whose stability has recently been 
analyzed El SOI El E] • While the existence of these compactly supported 
stable stationary states is a good indication of confinement properties for these 
repulsive-attractive potentials, it is not equivalent to confinement. Let us finally 
mention, that repulsive-attractive potentials have also been used in second order 
models for swarming |T8l HU [151 [12] where exponential decaying at infinity 
potentials are more suitable from the modelling viewpoint. 

Confinement properties were addressed in [lOJ taking advantage of the well- 
posedness theory for weak measure solutions of ([T]) developed in [TT]. Using the 
continuity with respect to initial data in the functional setting of probability 
measures ^(IR^), the authors reduced the confinement of the solutions to ([!]) 
to a similar confinement property for solutions of the associated particle ode 
system: 

Xj = — mj'VW{xi — Xj), i = 1, . . . ,n, (2) 

where Z{i) = {j e {!,..., : j ^ i,Xj{t) Xi{t)}, Xi{0) E for all 
i = l,...,n, < rrii < 1, and = M > 0. The authors obtain a 

confinement property in probability measures assuming that the potential is 
radial and attractive outside a ball, apart from other technical assumptions 
related to the well-posedness theory for probability measures. We will improve 
over the main result in [T^ in terms of the assumed attractive strength at 
infinity. More precisely, we will allow for slower growing at infinity potentials, 
see Section 2 for the precise hypotheses. 
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We will later obtain a confinement property for solutions to ([T]) in a smooth 
functional setting of compactly supported initial data in W^'°°(]R^). The trade- 
off is to allow more singular repulsive at the origin potentials. In fact, we will 
concentrate on the particular case of Newtonian repulsion plus smooth attractive 
potential with certain growth at infinity. In this functional setting, we can 
deal with smooth solutions obtained by a slight variation of the arguments in 
[HI in m E] . Section 3 shows that the same ideas used for particles and solutions 
in the functional setting V(M.^) apply in the continuum model ([T| for smooth 
solutions with this particular potential. The strong Newtonian repulsion at the 
origin of the potential allows us to derive a priori L°° bounds that otherwise are 
not known. Let us finally mention that confinement for the repulsive Newtonian 
plus an attractive harmonic potential was obtained in [9]. 

In the last section of the paper. Section 4, we perform some numerical com- 
putations for ^ with different repulsive-attractive potentials and we study the 
confinement of the stationary states for these cases. The stationary states are 
found using numerical techniques similar to the ones used in [HI EOj. We do 
a careful study in each case by computing the radius of the support in order 
to verify numerically the confinement for the solutions. Our numerical studies 
indicate that the assumptions we impose on the potentials are not sharp and 
could possibly be improved. 



2 Confinement for probability measures 

In this section, we will work with the theory developed in [11] in the framework 
of optimal mass transportation theory applied to ([T|. We remind the reader 
that the equation ([T]) can be classically understood as the gradient flow of the 
interaction potential energy [T3l [3l [H]. Optimal transport techniques allow 
to construct a well-posedness theory in the space of probability measures with 
bounded second moments 7^2(1^^) at least for smooth potentials [3j. The regu- 
larity assumptions on the potential were relaxed in [[llj allowing for potentials 
attractive at the origin with a at most Lipschitz singularity there, i.e., allowing 
for local behaviors like W{x) ^ with 1 < a < 2. 

More precisely, we assume that the potential W{x), see [TT], satisfies 

(NLO) W e C(M^) n C^(M^\{0}), W{x) = W{-x), and W{0) = 0. 
(NLl) W is A-convex for a certain A G M, i.e. W{x) — is convex. 
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(NL2) There exists a constant C > such that 

W{z) < C{1 + l^n, for all z eR. 

to derive the well-posedness theory of gradient flow solutions to Q with initial 
data in P2(R^). 

Under this set of assumptions (NL0)-(NL2), we can derive from |JJL, The- 
orems 2.12 and 2.13] that the mean-field limit associated to the model in ([T]) 
holds. On one hand, this means that approximating the initial data by atomic 
measures, we can approximate generic solutions of Q by particular solutions 
corresponding to initial data composed by finite number of atoms (particle so- 
lutions). On the other hand, this also implies that the solution of Q given 
in [m Theorems 2.12 and 2.13] coincides with the atomic measure constructed 
by evolving the locations of the atoms through the ODE system pi). In other 
words, if one is interested in showing a confinement property for (flF, it suffices 
to prove the confinement property for the particle system solving ^2]) since the 
solutions of the particle system (|2| approximate accurately in finite time in- 
tervals the solutions of the partial differential equation ([T]). All these details 
are fully explained in [ini Section 3] allowing us to reduce directly to particle 
solutions. 

To show confinement, we need additional assumptions on W{x) as in [10]. 
Throughout this paper, we assume that W is radially symmetric, and attractive 
outside some ball, i.e., 

(NL-RAD) W is radial, i.e. W{x) = w{\x\), and there exists Ra > 
such that w'{r) > for r > Ra- 

It is pointed out in [TTl Remark 1.1] that (NLl) guarantees that the repulsive 
force cannot be too strong, more precisely, 

Cw:= sup \VW{x)\ (3) 

X<^B(0,Ra)\{0} 

is bounded above. Here, we take the convention Cw = in case Ra = 0. 

In order to prove confinement results, we need some other condition to ensure 
that the attractive strength does not decay too fast at infinity. In addition 
to (NL0)-(NL3) and (NL-RAD), we assume that W satisfies the following 
confinement condition: 

(NL-CONF) lim w'{r)r = +oo. 
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which is less restrictive than the assumption in [TO], namely lim w'{r)^/r = +00. 



Therefore, our goal in this section is to show that if the particles interact un- 
der a potential satisfying assumptions (NL0)-(NL3), (NL-RAD), and (NL- 
CONF), have total mass M = 1, center of mass at 0, and are initially confined 
in -8(0, -Ro); then they will be confined in some ball -8(0, R) for all times, where 
R is independent of the number of particles n but only depending on the ker- 
nel W and the initial support of the cloud of particles Rq. Note that the zero 
center of mass assumption is possible due to the translational invariance of ([T]) 
and ([2]). Note also that the solution to the particle system ^ in the sense of 
[im Remark 2.1] might lead to a finite number of collision times, in which the 
solution may lose its regularity. Hence when we study the evolution of some 
quantities in time, we only take the time derivative in the time intervals in 
which the solution is regular. 

The strategy to get confinement for particles is as follows: we need to control 
quantities that quantifies how much the distribution spreads in time. In [TOt 
Proposition 4.2] the argument was based in following the particle furthest away 
from the origin and use some energetic arguments to control the mass of the 
particles nearby pushing the furthest particle. Here we follow a different idea. 
We consider other moments of the particle system to control the spread of 
the distribution of particles in conjunction with the evolution of the furthest 
particle from the center of mass. More precisely, we couple the evolution of the 
third absolute moment of the particle system with the evolution of the furthest 
particle. 



2.1 Evolution of the third moment 

Let M3(t) denote the third absolute moment of the particle system (|2]), namely 

n 
i=l 

In this section our main goal is to estimate the time derivative of M^lt). As we 
discussed before, there might be a finite number of collision times in which 
M2[t) becomes non-differentiable. Nevertheless, since all the particles have 
finite velocity, M^lt) is Lipschitz continuous in time even during collision. In 
all the computation below, the time derivative of M3 is only taken in the time 
intervals where M3(t) is differentiable; and the continuity of M^it) ensures that 
the fundamental theorem of calculus still holds for M^lt). 
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Since \x\^ is a convex function on M^, we know that every pair of attracting 
particles would give a negative contribution to dM^/dt, whereas every repulsing 
pair gives a positive contribution. We can directly evaluate dM^/dt as follows: 

dM^{t) 



dt 



2 ( -^i 






1 OC Ob j 1 
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1 X 'I OC j 1 





3 ^ ^ Tni \xi I (ij, Xi) 

i 

= 3^ ^ mimj{-VW{xi - Xj),Xi\xi\) 
i jez{i) 

= ~ 2 5^ 5Z ^i^jw'il^i - Xj\)Tij, 
i jez{i) 

I) and where antisymmetry of \/W{x) is 

I i j I 

used. Elementary manipulations yield that Tjj can be rewritten as 

Tij = {\xi \ + \xj\)- 
which gives the following upper and lower bound for Tij-. 

2 ( I I ~^ I '^j I ) I '^j I — '^'^i — ( I I ~^ I I } I I ■ ( ^ ) 

Next we will find an upper bound for —w'{\xi — Xj\) in Q. Let us define the 
nearest particles set N{i) as the set of indexes of particles that are possibly 
repelling the i-th particle, more precisely, 

N{i) := |j G {!,..., n} : < \xj{t) - Xi{t)\ < 

Then ^ implies that —w'{\xi — Xj\) < Cw for all j G N{i). 

For j ^ N{i), (NL-RAD) gives that —w'{\xi — Xj\) < 0. A better bound can 
be obtained using (NL-CONF): note that for any fixed constant i^i > to be 
specified momentarily, there exists some Rk^ > "iRa, such that 

—w'(r) < for all r > Rr,- 

r 

Let us define the set of furthest particles F{i) as the set of indexes of particles 
whose distance to the i-th particle are larger than Rk^, namely 



F{i) := <{j G {!,...,«} : \xj{t) - Xi{t)\ > Rr, 
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The definitions of N{i) and F{i) are illustrated in Figure [l| Then the upper 
bound for —w'{\xi — Xj\) can be summarized as following: 



w 



w'{\Xi - Xj\) < < 



for j G N{i), 
for j G F{i), 
for J ^ N{i)UF{i). 



(6) 




Figure 1: Illustration of the sets N{i) and F{i). For the i-th particle, N{i) is 
defined as the set of indexes of particles in the red region, while F{i) is the set 
of indexes of particles in the blue region. 

By plugging ^ and ^ into ^ and setting Ki := lOCwRa, we obtain 
^^^^f^ < mjCwRai\xi\ + \xj\) - mjKi^{\xi\ + 

< lcwRaJ2MT: -5n). (7) 



with 



Tr '■= ^ rnj{\xi\ + \xj\) and := ^ mj{\xi\ + \xj\) . 
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Now we claim that 

T;<4T^ ii\x,\>RK,. (8) 

Its vahdity is one of the main reasons for imposing the requirement (NL- 
CONF). To prove the claim, recall that we assume the center of mass is at 
at t = without loss of generality. Due to the conservation of the center of 
mass, for any time t, we have Xj{t) satisfies ^ ■ rrijXj (t) = 0. Let ei G denote 
the unit vector pointing in the direction of Xj, then it follows immediately that 

n 

rajXj ■ Cj = 0. 

For \xi\ > Rki, we split the above sum into three parts, and get 

^^^^^ TY\/ jOC j ' G-i^ — ^^^^^ TY\/ jOC j ' ^^^^^ TY\/jOCj ' 

3&N(i) j&F(i) j^N{i)UF{i) 

<-J2 ^3^1 ■ < Tl (9) 

where the first inequality is due to the fact that Xj ■ Cj > for all j ^ F{i). 

Moreover, recall that we find Rk^i we force it to be bigger than 2Ra. This is 
to guarantee that for all \xi\ > Rk^ > ^Ra and j G N{i), the angle between the 
vectors Xj and Xj is less than |. As a result, we have Xj ■ Cj > Noticing 
that for \xi\ > Rk^ and j G N{i) we also have \xj\ > \xi\/2, which is equivalent 
with \xj\ > |(|xj| + \xj\). Thus finally we have 

and by combining it with (|9]) we obtain the claim (|8|. 

Due to (|8|, we deduce that for any \xi\ > Rri, — 5T^ < — T*, hence ([T]) 
becomes 

with 

Ti := ^ Y mimj{\xi\ + \xj\) . 



Note that we can easily bound Ti by a constant only depending on W (since 
< Rki and \xj\ < Rk^ + thus we can rewrite the above inequality as 
inequality as 

where Ci,C2 only depends on W. At this point we will take a pause on the 
evolution of M3; we will revisit the inequality (10) soon in Section 2.2 to couple 
it with the evolution of the support. 



2.2 Coupling with the evolution of the support 

For all t > 0, let R{t) denote the distance of the furthest particle from the center 
of mass (which we assumed to be without loss of generality), namely 

R{t) := max \xi{t)\. (11) 

1=1,. ..,n 

It is pointed out in [111 Proposition 4.2] that R{t) is Lipschitz in time. Our goal 
is to prove that limsupi?(t) < R, where R only depends on W. 

t—>oo 

We begin by reminding a claim proved in [TUl, Proposition 2.2]: Let e be any 
unit vector. Then ^ 

J2 rnj>-J2mj, (12) 

Xj-e<R{t)/2 j 

i.e. the green region in Figure [2] contains at least 1/3 of the total mass. 

It is argued in the proof of Proposition 4.2] that for all time t > 0, there 
is a particle index io{t) (here zq may depend on t), such that 

r/+ T- (f) 

\x.M=m and ^^W=*^oW-^, (13) 

where ^ stands for the right derivative. This technical point is due again to 
the lack of regularity of R{t) for all t, see [HI Proposition 4.2]. From now on, 
the index io refers to any index satisfying the previous properties. 

To control ^R(t), it suffices to look at the outward velocity of the io-th 
particle at this time. As argued in [TOj Proposition 2.2] and illustrated in Figure 
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Figure 2: For any unit vector e, the green region above contains at least one 



third of the total mass. Here R{t) is as defined in (11) 



|3} the red region is possibly pushing it out, but all the green region is pulling it 
towards the origin. We proceed by estimating the compensation between these 
two competing effects. Let G{io) denote the set of indexes of particles in the 
green region in Figure [3} namely 

G(zo):= {jG{l,...,n}:x,(t)-^<^|. 

I 2^40 I 

Recall that throughout this section, we assume the total mass is 1 without loss 
of generality. It then follows from (|T2|) that 



E ^.>V3. (14) 

ieG(io) 

Using (|13|), g, and some simple manipulations, the growth of R{t) is con- 
trolled by the following inequality: 

—R{t)< Cw E rrij - E mjw'{\xj - Xig\)cos{9{-XiQ,Xj - Xig)), 
jeN{io) j&G{io) 

< Cw Yl ^^~\ nijw'{\x^-Xi^\), (15) 

jeN{io) ieG(io) 
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Figure 3: When R{t) > 2Ra, the particles in the green region are all pulling the 
particle Xi^ towards the origin, while the particles in the red region are possibly 
pushing it out. 

where 9{—Xif^,Xj — Xig) denotes the angle between the two vectors —Xi^ and 
Xj — Xig, which is less than | as shown in Figure 3j Due to (NL-CONF), for 
any large constant K2, which we will fix later, there exists some radius Rk2j 
such that w'{r) > K2/r for all r > Rk2- Hence for all t satisfying R{t) > Rk2, 



(15) becomes 



< Cw rrij - 



i6G{io) 



(16) 



12R(t) 



where (14) was used to obtain the last inequality. 



To ensure the coupling between the growth of M^{t) and the growth of R{t) 



go smoothly, let us go back to ( 10 ) and perform some elementary manipulation 
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on it. When R(t) > Rk^ + Ra, for any j G N{io), we have \xj\ > Rri, hence 



dt 



< C1-C2 ^ (^mi\xi\ ^ 



(17) 



And if in addition we have R{t) > 2(i?Xi+-Ra), then it fohows that G{io) C F{i) 
for any i G A^(io); hence we can replace the F{i) in (17) by G{io) and obtain 



ieAf{io) 



ieAf(io) 



rrii 



ieN{io) 



where we used (14) again to obtain the second inequahty. Finally, we set -Ri : = 
max{i?i^2, 2{Rki + Ra)}, to ensure that both (16) and (18) hold for R(t) > Ri. 

Finally we are ready to couple Msit) with R{t). By putting the estimates on 
^R{t) and ^M{t) together, we will show that if R{t) grows from Ri to some 
very large number in some time interval [ti,t2]; then the integral of dM^/dt 
is negative over this time interval, i.e., M3(t2) < M^lti). On the other hand, 
we will directly prove that M3(t2) must be bigger than M^iti), which causes a 
contradiction. 

Let Ai be a sufficiently large constant which we will determine later. If the 
particles start in B{0,Rq) and eventually touch the boundary of B{0, AiRi), 
then there exist < ti < ^2? such that 



R{h) = Ri, i?(t2) = -^Rit2)>0 



and 



R{t) G [Ri,AiRi] for all ti < t < ^2- 

More precisely, by letting ^2 := min{t > : R{t) = AiRi} > 0, and ti :- 
max{0 < t < t2 : R(t) = Ri}, they would satisfy all the requirements. 



Since R{t2) > R(ti), we have 

^d+ 



dt 



R{t))R{t) dt 



t2 



> 0. 
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Using (16), the above inequality implies 

ft2 



Cw I XI 1 ^(^) > /" dt, 

N{io{t)) J 



and by plugging it into the integral version of (|18|) we obtain 

(3 



*2 dM^ 



dt 



-(t) dt < 



dt. 



hence by choosing K2 := 1 + 36Cw{Ci + \C2Ra) /C2, which only depends on 
W, we have M^{t2) < M^^ti). 

On the other hand, if R(t) successfully grows from Ri to AiRi, we will show 
that M3 indeed has to increase, namely M^{t2) > 2M3(ti) for Ai sufficiently 
large. First, we can bound M3(ti) above by the very rough bound R(ti)^ = Rf. 



At time t2, recall that ^R(t2) > 0, hence (16) implies that 



^ ^ K2 



, , VlCvj A\R\ 



Finally, by noticing that 

12CyyAiRi 



Ms{t2) > (i?(t2) - Ra? Yl ^ 19r^4 fi> - ' 



ieAf(io(t2)) 

we obtain M3(t2) ^ (^i-Ri)^- Therefore we can choose Ai sufficiently large such 
that M3(t2) > 2M3(ti), which leads to a contradiction with M3(t2) < M3(ti). 

Note that the proof above shows that R{t) can never reaches R := AiRi, 
which is a large constant only depend on W, and in particular is independent 
of the number of particles. 



3 Confinement for kernel with Newtonian re- 
pulsion 

In this section, we consider the interaction kernel W{x) given by 

W{x) = -^{x) + Wa{x) (19) 
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with N >2. Here N'{x) is the Newtonian kernel, namely 

;^ln|a;| = 2, 
> 3, 



U{x) = {^^ Cn 



\x 



N-2 



where cn denotes the volume of a unit ball in M.^ . Throughout this section we 
assume that Wa{x) satisfies the following assumptions: 

(Wl) AWa e LL(M^). 

(W2) AWa is bounded in \ 5,(0) for any e > 0. 

(W-RAD) Wa{x) = w{\x\) with w G C\{0, oo)) and w'{r) > for r > 0. 

(W-CONF) lim w;'(r)r^/^ = +oo. 



Our goal in this section is to show that under the above assumptions, if a 
solution has total mass M = 1, center of mass at 0, and are initially confined in 
B{0, Ro), then it will be confined in some fixed ball centered at for all times, 
where the radius of the ball only depends on Wa, Rq, the dimension A^, and the 
norm of the initial data. 

Remark 1. For N = 1, the Newtonian kernel becomes \x\. Note that in this 
case the confinement result does not hold under the assumptions above, since the 
repulsive velocity field between two particles will be a constant regardless of the 
distance between them, while the attraction may vanish as the distance goes to 
infinity. We can compensate this difficulty by imposing stronger assumption on 
the attractiveness of Wa at infinity. More precisely, by replacing (W-CONF) 
by lim {w'{r) — l)r = +oo, the confinement result will hold with a similar proof 

r—^oo 

as in Section 2 carried over at the continuum level. 

Remark 2. (W-CONF) is more restrictive than (NL-CONF), especially for 
large N. In the proof below, one can see that dM^/dt does not cause a problem 
at all, indeed it satisfies the same inequality as in the non-singular kernel case 
in Section 2. The problem lies in dR/dt: due to the singular repulsive kernel. 



we got a worse control of dR/dt than before, see (25) 



We point out that slight variations of the arguments in [5, Section 5] and 
give a well-posedness theory for smooth solutions constructed by characteristics. 
More precisely, for any compactly supported initial data po ^ W^'°°(]R^), there 
exists a unique classical solution p G C\[0,T] x M^) n W/j,"^(M+, Wi'°^(M^)) 
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to Q with W satisfying (W1)-(W2). Moreover, the associated velocity field 
v(t,x) = —VW * p is Lipschitz continuous in both space and time, hence the 
characteristics are well defined: 

j^Xt = -{VW*p){t,Xt), 

and the solution p is given by 

pit,x) = po{X-')detiDX-'). 

Since the initial data is compactly supported, it remains compactly supported 
for all time (although the support may grow in time), and its support is obtained 
through the C^-characteristic maps Xt. 

First we remind a lemma showing L°°-bounds of the solution. This is a 
classical argument that can be seen for instance in [211 E] and O Section 5] but 
we give a short proof for completeness. 



Lemma 1. Let W be given by (19), with Wa satisfying (W1)-(W2). Let p be 
a classical solution to ([T]) with compactly supported initial data po e W^'°°(]R^). 
Then ||p(t, ■)||oo < Mq for allt>0, where Mq only depends on Wa and Pq. 

Proof. Due to the assumption (Wl), we can find tq > sufficiently small, such 
that 

/ \AWa{x)\dx<l. 
JB{0,ro) ^ 

Then it follows from (W2) that Mw := sup \AWa{x)\ is finite. We 
define Mq as 

Mo := max{2Mvy||po||i, IIPoIIoo}- 

Let us denote by M(t) = max^gjgjv p(t, x). If the desired result does not hold, 
then there exists some t\ > 0, such that M{ti) > Mq and M{t) is increasing at 
t = ti. This enables us to find some xi G M^, such that M{ti) = p{ti,xi), and 
^p(ti,xi) > 0. On the other hand, since p is a classical solution, we have 

d 

— p(ti, xi) = Vp ■ (WW * p) + piAW * p) 

dt (20) 

= p(tl, Xi) (^(p * AWa){ti,Xi) - p(ti, Xi) j . 



15 



Now let us split the integral p * APVa in 5(0, ro) and outside to get 

(p*Aiy„)(ti,xi) < ^||p(ti,-)||oo + M^||p(-,0)||i = ^p(ti,xi) + M^||poi|i. 



Plugging it into (|20|), we have 

1 \ /Mn n(t..T.A\ 

<0, 



d_ 
Ft 



p(ti,a;i) < p(ti,xi)(^MH/||po||i - 2^1^!' ^1) ) ^ PV'\^^\)\^ ^ 



which contradicts with the assumption that -^p{ti,xi) > 0. □ 

Next we present a technical lemma which will be used in the proof of confine- 
ment. 

Lemma 2. Assume u G L°°(M^) with < p < N . Then it follows that 

/ ^^^/y<C{N,p)\\u\\t!' [ u{y)dy] (21) 

Jb{xo,r) Fo — 2/r \Jb{xo,R) J 

for all xo G and all R> 0. 

Proof. First note that it suffices to prove the following inequality holds for all 



rfy<C(iV,p)||t;ri^||t;||i^-^)/^, 



by letting v{y) = XB{o,R){y)u{y + xq) , where xn is the indicator function on a set 
Q C M^. We point out that one could use Holder inequality and interpolation 
inequality on weak spaces to obtain a slightly weaker inequality than above, 
but we will use an easier and more elementary approach instead. 

Let w be an indicator function taking value ||f ||oo on some disk centered at 
and taking value outside, where the size of the disk is chosen such that w 
and V have the same norm. More precisely, w is given by 



w 



|^^||ooXB{o,ro), where Tq := 



^ 1 



Cn\\V\\oc 



l/N 



here ctv is the volume of the unit ball in M^. Since \\v\\i = \\w\\i, it is straight- 
forward to verify that 

/ v{y)dy < / w{y)dy for all r > 0. (22) 
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Now we start with the left hand side of (21), and Fubini's theorem yields that 

viy) 



Rjv \y\ 



-dy 



oo 



< 



JB(0,t~i/P) 
^Civ||f||oo| lAT. 

N — p 



^iy)X{t<\y\-p}dtdy 
w{y)dydt = 



JB{0,t~^/p) 

w{y) 



v{y)dydt 



m \yr 



-p Ljy II t'lloo 



■dy 



N -p 



\{N~p)/N 
ll 



where (22) was used. 



□ 



3.1 Evolution of the third absolute moment 

Similar to the particle system case, we also start with estimating the time 
derivative of the third absolute moment M3. Here the third absolute moment 
M3 is given by 

Ms^t) := / p{t,x)\x\^dx, 



and note that in the continuum setting M3 is indeed differentiable in time for 
all t >0, since p{t, x) is a classical solution. The same computation as ^ leads 
to 

dM^it) 3 f f , x — y I IN , , 

- p{t,x)p{t,y)w {\x -y\)- : ■ {x\x\ - y\y\)dydx. 



dt 2 J^N J^N \x — y\ 

Due to the assumptions (W-RAD) and (W-CONF) on Wa, for any A > 
(which will be fixed at the end of this subsection), there exists some Ra > 1, 
such that the following bound for —w'{r) holds, where Cn is some constant only 
depending on N: 



-w'{r) < 



'^^ for < r < Ra, 



for r > Ra- 



-l/N 
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Using this bound and ([s]), ^^^^^ becomes 



dMsjt) ^ 3 
dt - 2 



b{x,Ra) F ^ y 



I 



x,Ra) 



p(t,y)^\x-y\^''~'^^''i\x\ + \y\)dy]dx. 



'Y'x 



Similar to Section 2, we again claim that < T^/2 for |x| > 2Ra. We start 
with controlling T^. It is easy to check that 

Jb{x,Ra) F ~ 1/1 

Note that the singularity of the Newtonian kernel is more difficult to treat than 
in Section 2. We compensate this difficulty by using the fact that pit,x) is 
uniformly bounded by Mq from Lemma [T} Hence for > 2, we are able to 
apply Lemma [2] to p(t, ■), and obtain 

T; < Cs{2\x\+Ra) ( [ p{t,y)dy) ' (23) 

\JB{x,Ra) J 

here C3 only depends on and Mq as obtained in Lemma [1} 

To simplify notation, from now on, we define by m(t, the mass of p within 
radius Ra of 2; at time t, namely 



m{t,x) := / p{t,y)dy. 

JB{x,Ra) 

Recall that in the beginning of this section we assume that po integrates to 1, 
which implies that p(t, ■) also integrates to 1 for all t > 0. Then, for any t 
and X, one of the two following scenarios must be true: either m{t,x) > |, or 
1 — m(t, x) > ^. 

If m(t,x) > ^ at some |x| > 2Ra, it follows that m(t,a;)^/^ is comparable to 
m(t,x). Hence, using (23), we get 

< 6-32^-2/^ m{t,x){2\x\+RA) 

< 8C3 / Pit,y){\x\ + \y\)dy (since |x| > 2Ra) . 

Jb{x,Ra) 
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Hence by repeating the same argument on the center of mass as in Section 2 
(see ^ and the paragraph after it), we can choose A to be sufficiently large, 
then we would obtain that < 

On the other hand, if the opposite scenario is true at some |x| > 2Ra, i.e. 

1 



/ 



p{t,y)dy > 

B(x,Ra) 



then one can directly bound by C(Vr, A^)|x| by applying Lemma [T] and 
Lemma |2} Meanwhile it follows directly from the definition of that > 
hence by choosing A sufficiently large we obtain that < Ta/2. 



Finally, we choose A to be the maximum value needed in the two scenarios. 
As a result, < Ta/2 holds for for all |x| > IRa-, implying that 



dt 



<Ci-C5 [ P{t,x)\x\ [ p{t,y)dydx, (24) 

jR^\BiO,2RA) Jr^\B{x,Ra) 



where 6*4,(75 only depends on Wa, N and ||po||oo- Note that this inequality is 
parallel to the inequality (10) for the discrete case. 



3.2 Coupling with the evolution of the support 



Next we will proceed similarly as in Section |2.2[ where most of the arguments 
are parallel. We will quickly go through the similar parts in the proof, and 
emphasize the differences caused by the Newtonian repulsive kernel. 

At time t, we can find Xq G 9supp(po) depending on t, such that 

|X,(a;o)| = R{t) and ^ < -{VW * p){Xt{x,),t) ■ , 
similarly to [3 H] . 

Due to (W-CONF), for any large constant to be determined later, there 
exists some radius Rk^ > 6Ra such that w'{r) > K^/r^/^ for all r > Rr^.- 
Hence whenever R{t) > Rk^, the growth of R{t) is now controlled by 

dR{t) ^ f CNp{t,y) , K: 



< 



dt - jBix,(x,),RA)\y-xM^W-' i2R{ty/^ 
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where the second term on the right hand side is obtained in the same way 
as the last term in (16), except that the power 1 is replaced by 1/A^ due to 
(W-CONF). 

To deal with the singularity in the first term on the right hand side, recall 
that ||p(t, OIloo is bounded above by Mq for all time due to Lemma [l| which 
again enables us to apply Lemma [2] to obtain 



dR{t) 
dt 



<C6m(t,Xi(xo))^/^- 



i2R{ty/^' 



(26) 



where Cg only depends on A^, Wa and ||po| 



Similar to Section 2.2, we can find some time interval [^1,^2], such that R{t) 
increases from Rk^ to A2RK3 within [^1,^2], and R{t2) > 0. Here A2 is a 
sufficiently large number to be determined at the end of this subsection. Then 
we have 

^R{ty/^ > 0, 

dt ^ ' 



t2 



implying that 



t2 



C,m{t,Xt{x,)Y"'R{tf"' 



K3 
12 



dt > 0. 



(27) 



We apply Holder's inequality on (27), and obtain that 

i2 / \ N 



m 



{t,Xt{xo))R{t)dt > 



12a 



(t2-tl). 



(2^ 



Note that this extra step is needed here but unnecessary in Section 2, due to 
the different powers in (NL-CONF) and (W-CONF). 



Now we are ready to couple the growth of M3 with (28). Since R{t) > 6Ra for 
all t G [^1,^2] (recall that when defining Rk:^ we set it to be greater than QRa), 
we could treat (24) in the same way as we did in (17) and (18), and bound the 
growth of M3 as follows: 

dMsit) 1 



dt 



< iC^ + -C^Ra) - C,m{t,Xtixo))Rit). 



(29) 



Then we integrate (29) in [^1,^2]; and it becomes 



"42 



C,m{t,Xt{xo))R{t) - (C4 + -C,RA))dt < M^{t{) - M^it^). 
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By putting the above inequality together with (28), we can fix to be suffi- 
ciently large such that M^{ti) > M3(t2)- 



Finally, we prove that if A2 is sufficiently large, we would have 71^3(^2) > 



M^(ti), hence causing a contradiction. It follows from (26) and R{t2) > that 



implying that 



M3(t2) > m{t2,Xt,{xo)) {A2RK, - RksY > AlRl 



X3' 



which can be made to be greater than Ms{ti) if A2 is chosen to be sufficiently 
large, thus we obtain a contradiction with M^lti) > M^{t2)- This means that 
R{t) can never reach A2Rk3-, thus implies the confinement of support for all 
times. 



Remark 3. Let us emphasize that for potentials W given by (19), we are only 
able to prove confinement in the continuum setting, not in the particle setting. 
The reason is that in the coupling method we use, we need to bound the repulsion 
part of ^R{t) using the mass in some neighborhood of the outermost particle. In 
the continuum setting this is achieved by first obtaining an bound on p{t, •) 
in Lemma^ then applying Lemma^to arrive to (26). However, in the particle 
setting, we are unable to obtain a bound on the "local density" of the particles 
that is independent of the particle number, and we are unaware of any such 
results for repulsive- attractive kernels to the best of our knowledge. Intuitively 
we do expect the "density" of particles to be bounded, since the singular repulsion 
would not allow the particles to be densely concentrated. We find it an interesting 
open problem to prove some non-local version of Lemma^for the particle system 
(0) with W given by (|l9[ . 



4 Numerics 



In this section we numerically check the confinement properties of several poten- 
tials together with the long-time behavior of the corresponding particle systems. 
Let us remark that in all the cases we have simulated, for which confinement 
holds, the long time behavior of the system seems to converge toward a com- 
pactly supported stationary state. In some of the potentials below, this has not 
been rigorously proved. This is an interesting theoretical question that will be 
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treated elsewhere. Our objective in this section is to check if the conditions 
under which confinement has been shown in previous sections are sharp or not. 
With this aim, we remind, as it was said in Section 2, that equation Q is a 
gradient fiow of the interaction energy 

E[p\ = \ I I W{x - y) p{x) p{y) dx dy 

with respect to the Wasserstein distance. Thus, stable stationary states of ([T]) 
are local minimizers of the interaction energy. 



In Section 2 we have shown that the radius -R(t) defined by (11 ) is bounded by 
a constant R that depends only on the potential W and the initial data. More- 
over, R is independent of the number of particles and under certain additional 
assumptions, see Section 2, we know that the particle systems are indeed good 
approximations of the solutions to the continuum model ([T]). For this reason, 
we have chosen a particle framework to perform our numerical investigation. 
We also follow the idea of decreasing the energy since stationary states are local 
minimizers of the energy. Given n particles located at xi, . . . , a;„ G with 
masses m\ = m2 = ■ ■ ■ = rrin = 1/n, their discrete interaction energy is given 

by 

^ n n 

E[xi, • • • , = ^ ^ ^ W{xi - xj). 
i=i j=i 



The simulations are done by an explicit Euler scheme leading to a trivial 
gradient descent method as long as the energy is decreasing at each time step. 
This method allows to efficiently solve for stationary states of ([2]). In stiffer 
situations, as for the Morse potentials below, an explicit Runge-Kutta method 
is used instead. These methods are essentially the same as the ones used in 
[1T| HO] for finding stationary states of different repulsive-attractive potentials. 
Our stopping criterion is to achieve a numerical steady state. For us, a numerical 
steady state is a particle distribution for which the discrete /°°-norm of the 
velocity field in ^ is below some predetermined threshold, which we impose to 
be 0.001/n. 

The section is divided into three subsections, each one showing the results 
for a particular chosen potential. The limit growth for the attractiveness of 
the potential at infinity under condition (NL-CONF) is log(r). For this rea- 
son in Section 4.1 we constructed a piecewise potential with exact logarithmic 
attraction at infinity. This selection has been done to check the sharpness of 
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condition (NL-CONF). In Subsection 4.2 we go further to take a piecewise 
potential growing at infinity exactly like log(log(r)), which grows even slower at 
infinity compared to log(r). This potential does not satisfy the condition (NL- 
CONF). At the end of this section, in Subsection 4.3, we also analyze the case of 
the Morse potential. This potential is known to be a repulsive-attractive poten- 
tial under certain choices of the parameters with negligible attractive strength 
at infinity, i.e., W{x) — )■ as — )■ oo. These potentials are more interesting in 
terms of biological relevance as discussed in [THl [12] • 

As a final remark, we point out that all the used potentials are not singular 
at the origin and simulations are performed in dimension N = 2. 



4.1 Logarithmic attraction at infinity 

We show confinement when the potential has exact logarithmic attraction at 
infinity. The chosen potential is 



wir 



95 2 _ 
2 ' 

log (r) 



83 



64r3 + 2|9^4 _ ^^5 0<r<l, 

r > 1. 



(30) 



It can be checked that it is a repulsive-attractive satisfying w{0) = w{l) = 0, 
w{r) G C^(0,-|-oo), and the repulsion at the origin is ~ — r. The stationary 
states are shown in Table [1} We have chosen initial data in such a way that 



n = 100 


n = 200 


n = 500 


n = 1000 










r ~ 0.2286 


r ~ 0.2315 


r ~ 0.2343 


r ~ 0.2352 



Table 1: Stationary states and radius of their support as a function of the 



number of particles for the potential w{r) given in equation (30) 



the particles feel the logarithmic interaction by randomly placing n particles in 



a centered square in such a way that \xi 



> 1 for some values of i,j e 



{1, 2, . . . , n}. We have run simulations varying the number of particles n and 
the initial data. For each simulation the center of mass Cn for the particle 
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system was computed and then 

r„(t) = max Ix.ft) — CA- 

It is observed that r„(t) < 1 in all the cases for large times and it converges to 
some asymptotic value. This fact can be explained because when all the particles 
are out of the range of the log(r) part then the radius and the behavior depends 
only on the polynomial part of the potential. Simulations indicate that there is 
confinement for this potential and thus, condition (NL-CONF) is not sharp. 
Figure |4] shows the evolution of the radius as a function of the particle number 
n and as a function of time for a particular initial data. 




(a) (b) 
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0.2 




0.1 




n 





200 400 600 800 1000 

n 



(c) 

Figure 4: (a) Evolution of the radius as a function of t for n = 1000. (b) Zoom 
of picture (a), (c) Evolution of the radius as a function of n. 
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4.2 Log-log attraction at infinity 



In this case we consider a potential behaving hke log(log(r)) at infinity 

+ 



wir] 



^ _ 2 r(r-e)^ ^ Kr-ef Q < r < e, 



log (log (r)) 



r > e. 



(31) 



This potential satisfies w{0) = w{e) = 0, it is C^(0, +oo), and the repulsion at 
the origin is ~ — r. The numerical experiments suggest that by increasing the 
number of particles, the radius of the support of the stationary state increases 
and stabilizes. Table [2] shows the stationary states as a function of the number 
of particles n for this potential. These numerical simulations, together with the 
evolution of the radius of the support both in time and as a function of the 
number of particles not shown here, indicate that even if the growth at infinity 
of the potential is less than log(r) there is still confinement. 



n = 100 


n = 200 


n = 500 


n = 1000 










r ~ 0.7838 


r ~ 0.7954 


r ~ 0.8052 


r ~ 0.8085 



Table 2: Stationary states and radius of their support as a function of the 



number of particles for the potential w{r) given in equation (31) 



4.3 Morse potential 

The usual form of this potential is the following 

t/(r) = -C^e-'^/^-^+C^e-^/^«, 

where constants Ca and Cr are the attraction and repulsion strength respec- 
tively and the constants Ia and Ir are their respective length scales. For our 
simulations we will take the scaling shown in [181 Il2]- That is, 

U{r)=CAiVir)-CV{r/l)), 

where V{r) = — exp(— r//^) and C = Cr/Ca and / = Ir/Ia- It is known for 
this potential [12] that for C > 1 and / < 1 the potential U{r) is short- 
range repulsive and long-range attractive with a unique minimum defining a 
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typical distance between particles. Also, in this regime, the condition Cl^ = 1 
distinguishes between the so-called H-stable and catastrophic regimes. 

H-stable case: In our simulations we fix the parameters as Ca = = 1? 
Cr = 1.9, and Ir = 0.8 leading to C = 1.9, / = 0.8, and Cl^ = 1.216 > 1. 



n = 100 


n = 500 


n = 1000 


n = 2000 










r ~ 15.96 


r ~ 34.002 


r ~ 47.35 


r ~ 64.45 



Table 3: Stationary states and radius of their support as a function of the 



number of particles for the potential U{r) given in equation (31). 



Numerical experiments in Table [3] demonstrate that the radius increases by 
increasing the number of particles, but with a slower rate. As clearly visualized 
in Figure |5](a), the radius appears to grow hke a square root function as the 
number of particles increases. This observation is further supported by numer- 
ical evidence in Figure [5|^b), where we plot the square of the radius versus the 
number of particles, and the linear regression provides a good fit to the data. 

It would be interesting to study the H-stable case in more details, although 
it is outside of the scope of this paper. From these numerical results, we can 
extract a conjecture that for the continuum system the support of the density 
would go unbounded over time, i.e., the confinement result should not hold for 
the H-stable potential. 

Catastrophic case: The parameters we choose for the experiments are Ca = 
= I^Cr = 1.3 and Ir = 0.2 so that CP = 0.052 < 1. 

The results for this case are shown in Table HI In contrast to the H-stable 
case, the radius of the support converges to a limiting value. In Figure |6]^a) we 
observe how the radius of the support decreases in time to a limiting value with 
n = 1000 particles, and in Figure |6](b) we show how the radius of the support 
of the stationary state increases and converges to a certain value as a function 
of the number of particles. We conclude that there should be confinement 
properties for the Morse potential in the catastrophic case. 

The final goal would be to find replacements for the condition (NL-CONF) 
in order to include the cases where w{r) — as r — t- oo. One possibility is 
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(a) 



(b) 



Figure 5: H-stable Case: (a) Evolution of the radius as a function of n. (b) 
Squared radius of the support of the steady state as a function of n and the 
hnear regression curve y = 170.34 + 1.99n computed using the crossed points. 



n = 100 


n = 500 


n = 1000 


n = 2000 










r ~ 0.5269 


r ~ 0.5480 


r ~ 0.5518 


r ~ 0.5540 



Table 4: Stationary states and radius of their support as a function of the 



number of particles for the potential U{r) given in equation (31). 



to invoke scaling limits for integrable potentials. More precisely, we scale the 
potential in (jl]) as 



(32) 



in such a way that We{x) = e '^Wlx/e) approximates a Dirac Delta at with 
certain weight as e — )■ 0. Now, if the potential is such that 



W{x) dx . 



then equation (32) is formally approaching 

= aV ■ (pVp), X G M^, t > 0. 
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Figure 6: Catastrophic Case: (a) Evolution of the radius as a function of t for 
n = 1000. (b) Evolution of the radius as a function of n. 



In the H-stable case, the Morse potential satisfying a > leads to a limiting 
nonlinear diffusive equation, which is coherent with the no confinement prop- 
erty. In the catastrophic case, the Morse potential satisfying a < leads to a 
limiting anti-diffusive nonlinear equation, which might also be coherent with the 
confinement property of the potential. We conjecture these integrability condi- 
tions might have some implications for confinement properties of potentials. 
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